Nanoscale heat flux between anisotropic uniaxial media 
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We present a theoretical study of near-field heat transfer between two anisotropic materials sep- 
arated by a small vacuum gap and maintained in a stationary non-equilibrium thermal situation. 
By combining standard stochastic electrodynamics and the Maxwell- Garnett description for effec- 
tive media, we show that heat flux can be significantly enhanced by air inclusions. This result is 
explained by : (a) the presence of extraordinary surface waves that give rise to supplementary chan- 
nels for heat transfer throughout the gap, (b) an increase in the contribution given by the ordinary 
surface waves at resonance, (c) and the appearance of frustrated modes over a broad spectral range. 

PACS numbers: 44.40.+a, 78.20.-e,03.50.De 



I. INTRODUCTION 



present our conclusions. 



Near field heat transfer—^ between closely spaced 
isotropic media has been intensively studied since it has 
been predicted that the heat flux at nanoscale can exceed 
the far-field limit of the Planck's blackbody theory by or- 
ders of magnitude 5,6 . When considering dielectrics, sur- 
face phonon polaritons provide additional enhancement 
as discussed in Refsi"— . Several experiments have re- 
cently confirmed the theoretical predictions for simple 
systemsiiii^"— . 

With the modern techniques of nanofabrication it is 
now possible to explore a whole new level of complexity 
in material science and to fabricate artificial materials 
that can exhibit a considerable diversity of optical prop- 
erties^—. In many situations, these composite media 
possess privileged orientations so that their electromag- 
netic response depends on the direction of photons propa- 
gation. When the photon's wavelength in such a medium 
is large compared to the size of its representative unit cell, 
the latter behaves effectively like an anisotropic material 
and therefore may be described by an effective permit- 
tivity tensor (and, when necessary, an effective perme- 
ability as well). This naturally points to the question of 
how anisotropy influences the near-field heat transfer. 

In this work, we address this question in the particu- 
lar case of two semi-infinite uniaxial media characterized 
by optical axes orthogonally oriented with respect to the 
surface of interaction. The paper is organized as fol- 
lows: In Sec. II we derive the expression for the heat flux 
between two anisotropic media. After a brief descrip- 
tion of the relevant composite media to our purposes in 
Sec. Ill we investigate in Sec. IV the surface and Brew- 
ster modes supported by them and their main features. 
Next, we compare in Sec. V the near-field heat exchanges 
between two uniaxial media to the classical ones between 
two isotropic media. Finally, in order to explain the dif- 
ference in the behavior of isotropic and anisotropic ma- 
terials, we discuss in Sec. VI the transmission factor in 
detail between two uniaxial media and in Sec. VII we 



II. RADIATIVE HEAT TRANSFER BETWEEN 
ANISOTROPIC MEDIA 

Let Bi and B2 be two anisotropic semi-infinite bod- 
ies, filling respectively the regions z < and z > d and 
leaving a vacuum gap between them. In order to ensure 
a stationary process, we assume that the Bi are in local 
thermal equilibrium at a temperature T^, with T\ 7^ T2. 
The heat flux P between the two bodies is given by 



P(Ti,r 2 ,a)= / dA-(S), 

JA 12 



(1) 



where S = E x H is the Poynting vector and A\2 is any 
surface that separates the two bodies. By taking such 
a surface to be a plane defined by z = zo (0 < zo < d) 
and using the (transverse) translational invariance of our 
system, the previous equation simplifies to 



P(T 1 ,T 2 ,a) = A{S z ), 



(2) 



showing that only the z-component of the Poynting vec- 
tor is needed. After a straightforward calculation, the 
latter can be conveniently written as^ 



{s. 



(3) 



where we identify the mean energy of a harmonic oscil- 
lator 



e(w,r) 



(4) 



and also the averaged spectral Poynting vector— 
(S u ) = 2ReTr J dv\ ^G(r,r / )0 z ^G t (r,r') 

-a z G t (r,r')^G(r,r') 



(5) 
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where r = rn + zz and G(r, r') is the electrical Green's 
dyadic, satisfying 



-e(r,cj) 



G(r,r / ,w)=(5(r-r / )l (6) 



Moreover, we have introduced Boltzmann's constant /cb, 
Planck's constant 2irh; the f symbolizes hermitian con- 
jugation and Tr the 3x3 trace. 

In order to evaluate the heat flux in the given geom- 
etry we have to determine the Green's dyadic inside the 
gap region. This can be done by considering the multi- 
ple scattering of a plane wave due to a source inside the 
gap22. Details and the final expression for the Green's 
dyadic can be found in appendix O When inserting the 
resulting expression in Eq. (|A13|) into the heat flux for- 
mula, we find 



(S u ) 



I 



d 2 t 



(2*)- 



■T(lo,k; d). 



(7) 



The integral is carried out over all transverse wave vectors 
k = (k x , k y Y including propagating modes with n < uo/c 
and evanescent modes with n > uj/c, where c is the ve- 
locity of light in vacuum. The energy transmission coef- 
ficient T(cj, n; d) is different for propagating and evanes- 
cent modes and can be stated as 



T(u,K]d) = 

f Tr[(l - R^R 2 )E> 12 (1 - RlRi)E> 12t ' 
[Tr[(f4 - R 2 )© 12 (Ri - lEll)© 12 ^- 



, K < Uj/c 

2 ^ d , K>UJ/C 

(8) 



From Eqs (|3]), ([7]) and (|5} we see that, once the reflec- 
tion matrices are known, it is possible to determine the 
heat flux between two arbitrary anisotropic semi-infinite 
bodies kept at fixed temperatures T\ and T 2 . Moreover, 
in order to have an independent check, we verified that 
Eq. ([7]) also can be derived from the general scatter- 
ing formalism derived on Ref^. In the following we will 
use these expressions to discuss the heat flux between 
two uniaxial anisotropic materials with their optical axes 
normal to the interface. 




Figure 1: Sketch of two porous slabs with different tempera- 
tures separated by a vacuum gap. 



where 7 r = s/oo 2 /c 2 — k 2 and Hi , R 2 are the 2x2 reflec- 
tion matrices characterizing interfaces. By writing them 
a bit more explicitly, 



R ? ; 



7f' S (cJ,ft) 7f' P (cJ,ft) 



(9) 



we see that their elements r^' A are the reflection coeffi- 
cients for the scattering of an incoming A-polarized plane 
wave into an outgoing A'-polarized wave. In the isotropic 
limit they reduce to the usual Fresnel coefficients 



»?>,«) = 



7r ~ \/ei(Lj)L0 2 /c 2 - K 2 

7r + V 'ei(w)w 2 / 'c 2 - K 2 ' 



(H&Yh - y/ej(uj)uj 2 /c 2 - K 2 

£i(w)7r + \/ £i(u)u) 2 / C 2 - K 2 ' 



r?' p (w>«) =rr(w,K) =0 



(10) 



and we see that the matrices become diagonal. In addi- 
tion, we have introduced the matrix D 12 , defined by 



D 



12 



(1 - RiR 2 e 2i7rd ) _1 , 



(11) 



which gives rise to a Fabry- Perot- like denominator for 
T(cj, d) in the isotropic case. 



III. POROUS MATERIALS 

The structures investigated in this paper are depicted 
in Fig. [TJ They are two semi-infinite media composed by 
a host isotropic material, defined by its complex dielec- 
tric function e^uj) = e f h (u) + ie^cj) (where e^((j) > 0), 
with uniform cylindrical inclusions oriented in the direc- 
tion orthogonal to the surface as shown in Fig. [TJ These 
inclusions in turn are filled by a medium of dielectric per- 
mittivity ei, that is also assumed to be isotropic. When 
the size of the representative unit cell is much smaller 
than all the other characteristic scales involved, a suit- 
able volume average of the material's local electromag- 
netic response can be made. In our case, the emerging 
azimuthal symmetry in this long wavelength limit gives 
rise to effective uniaxial crystals with a permittivity ten- 
sor of the form 



e en 



[e x <g> e x + e y <g) e y ] + e±e z <g) e z 



(12) 



where e x , e y , and e z are orthogonal unit vectors in x, y, 
and z direction. The parallel and perpendicular compo- 
nents can be derived from the Maxwell- Gar net t effective 
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medium theory (EMT) 2 ^ 



25,26 



(a) 



£|| = £h 



ei (l + /) + Ch (l-/) 
£i (l-/) + e h (l + /)' 
e h (l - /) + ej, 



(13) 
(14) 



where / is the volume fraction of inclusions. For the 
structure considered in this work the deviation from the 
exact result of homogenization given in Refs j 28 i 29 is small 
even for relatively high filling factors such as / = 0.5. 
Hence, we will discuss the heat flux between porous me- 
dia with the Maxwell- Gar nett expression for / G [0,0.5] 
in this work. 

The condition of long wavelengths sets a limit to the 
lattice constant a of the inclusions for which the EMT 
can be used. In the far-field regime this condition is 
fullfilled when the thermal wavelength A t h = hc/k^T is 
larger than a. In the near-field region the contributing 
modes at a distance d above the porous material have 
a lateral wavelength which depends on d. For n = 2tt /a 
(which corresponds to a lateral wavelength a) the evanes- 
cent waves are damped as exp[— yj (2it / a) 2 — uj 2 /c 2 d] « 
exp(— (2k / a)d) in the non-retarded near-field region 
above the porous material. It follows that the contri- 
bution to the heat flux is dominated by evanescent waves 
with lateral wavelength larger than a if d > a/(2it). On 
the other hand, one can argue that a nonlocal model for 
the permittivity is necessary if the lateral wave vectors 
k are on the order of ir/a. Since the exponential in the 
transmission coefficient in Eq. ([8]) for k > uj/c sets a cut- 
off for Hi of the modes contributing to the near-field heat 
flux which is ~ one finds that a local EMT descrip- 
tion is permissible if d > cl/tt. Hence, for a given lattice 
constant a of the inclusions the validity of the EMT in 
Eq. ([T3]) in the near field regime is restricted to d > a/ir. 
Artificial structures as depicted in Fig. [T] can have an a 
on the order of 100 nm^ 2 . so that the distances for which 
the EMT can be considered as appropriate in this case 
are about d > 30 nm. Nonetheless, chemically produced 
nanoporous materials can show smaller structures 2 ^ so 
that we will consider distances d G [10 nm, 100 /im]. 



(b) 
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IV. SURFACE AND BREWSTER MODES IN 
POROUS MEDIA 



Figure 2: Plot of the dispersion curves (white dashed lines) 
from Eq. ([22]) in the (lj,k) plane for filling factors (a) / = 0.1, 
(b) / = 0.3, and (c) / = 0.5. The white dash-dotted line 
represents the light line in vacuum (uj = kc). Furthermore 
the dark (blue) areas mark the region for which 7 P is purely 
real, whereas the bright (red) areas are the regions for which 
7 P is purely imaginary. 



Let us study the surface waves supported by these me- 
dia when they are sufficiently far away from each other so 
that any coupling of evanescent waves can be neglected. 
By definition, these surface waves are resonant surface 
modes and therefore are determined by the poles of the 
reflection coefficients of these media. For out-of-plane 
uniaxial media the components of the reflection matrix 



r P, P (^^) 



s,p — ' p,s 



7r ~7s 
7r - 7s ' 
. € ||7r ~7p 
C||7r + 7p' 

r n * = 0, 



(15) 

(16) 
(17) 
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Figure 3: Plot of ln(l/|r p , p | 2 ) in (u),k) plane for (a) / = 0.1, 
(b) / = 0.3, and (c) / = 0.5. 



where 7 S5P are given by the solutions of Fresnel equations 
in the anisotropic material 30 



7s 



€\\UJ 2 /c 2 — ft 2 , 



7 P = ^e\\uj 2 /c 2 



(18) 
(19) 



and hence it follows at once that the surface modes are 
determined by 



(7r + 7s) = 0, 
0||7r + 7p) = 0. 



(20) 
(21) 




It is straightforward to verify that in this case only the 
second equation above can be satisfied, meaning that 
only p-polarized surface waves can exist at the interface 
of these media. Solving that equation explicitly for n 
gives us the sought dispersion relation of surface waves 



(22) 



but one must be aware that ([22]) has two branches, and 
only one is connected to surface modes 39 . Since their 
dispersion relation involves en and e_L, these waves are 
also called extraordinary surface waves 3 -^, and they reflect 
the material anisotropy. When en = e± = e, Eq. ([22]) 
degen erates int o the well-known dispersion relation n = 
uo/c^/e/{e + 1) of surface modes supported by a semi- 
infinite isotropic medium (bounded by vacuum) with a 
dielectric permittivity e. In Fig. [2] we plot the dispersion 
curves for silicon carbide (SiC) with vacuum inclusions 
for different filling factors / = 0.1, / = 0.3 and / = 0.5. 
The dielectric function of SiC is described 3 -^ by the simple 
model 



eh 



iouT 



(23) 



where uj l = 1.827 • lO^s" 1 , uj t = 1.495 • 10 14 s-\ 
T = 0.9 • 10 12 s _1 , and = 6.7 denote respectively 
the longitudinal and transversal optical phonon pulsa- 
tion, the damping factor and the high frequency dielec- 
tric constant, respectively. In order to avoid the inherent 
difficulties of multiple possible interpretations of complex 
dispersion relations 33 , we have deliberately neglected the 
host material losses to represent these curves. The rele- 
vance of this approximation can be checked by comparing 
Fig.[2]with Fig.[3j where we plot the reflection coefficients 
of dissipating porous material. In order to distinguish 
between evanescent and propagative waves inside the ef- 
fective medium, solutions of Eq. ([22]) are superimposed 
in Fig.[2]to a two-color background. This background is a 
binary representation of ( = sgn(e||o; 2 /c 2 — e\\K 2 / e . In 
the blue zones ( < so that only evanescent modes can 
exist, and conversely, in the red zones we have £ > and 
all modes are propagative. Similarly the light line uo = ck 
allows us to distinguish between the radiative (propaga- 
tive) and the non-radiative (evanescent) modes inside the 
vacuum. Notice that, in order to satisfy Eq. ([22]) both £ 
and sgn(cj 2 /c 2 — n 2 ) must be the same. In other words, 
frustrated modes cannot satisfy the dispersion relation 



Now let us turn to the description of modes supported 
by our artificial structures. For low filling factors we 
note in Fig. [2] (a) the existence of two surface modes. 
The first one (at a lower frequency) corresponds to the 
classical surface phonon-polariton (SPP) supported by a 
massive SiC sample 34 . That surface mode is also present 
in isotropic SiC. The most interesting feature of Fig. [2] 
(a) is, however, the appearance of a second surface mode 
at higher frequencies, because it is a signature of the 



5 



anisotropic character of the material and therefore a di- 
rect consequence of the vacuum inclusions in the host 
medium. As the porosity increases, both surface waves 
split. Beyond a critical filling factor between / = 0.3 
and / = 0.5, the upper surface wave disappears as is 
seen in Fig. 03 Nevertheless the SPP which still exists 
continue to move toward the smaller frequencies, i.e., to 
ujt- Above the light line, we see that the anisotropy 
gives rise to two different types of Brewster modes. At 
high frequency we recognize the usual modes where the 
reflection coefficient [Fig. [3] (a)] of the effective medium 
vanishes. In addition to these modes, different Brewster 
modes appear depending on the value of filling factor. 
Also, we see on the reflection curves (Fig. [3]) that the 
Christiansen point 35 for which the reflectivity is zero for 
all n does not depend on the porosity. Indeed, an inspec- 
tion of expressions (p~3|) and ([T4]) shows that the condition 
for the Christiansen point of the host material eh = 1 im- 
plies that e\\ = 1 and e±_ = 1 so that, according to (fT7|h 
the reflection coefficients vanish. 



V. HEAT FLUX BETWEEN POROUS MEDIA 

Before we discuss the influence of the inclusions on 
the heat flux, we show in Fig. |4] the results of the mean 
Poynting vector (S z ) between two semi-infinite SiC bod- 
ies at fixed temperatures Ti = 300 K and T2 = K. First 
of all one can see that the heat flux becomes very large 
for distances much smaller than the thermal wavelength 
A t h = hc/k^T (which is about 7.68 /im for T = 300 K). 
At d = 10 nm the heat flux for the two SiC bodies is 
about 1000 times larger than the heat flux between two 
black bodies. This increase is due to the frustrated to- 
tal internal reflection and to the coupled surface phonon 
polariton modes 9 . In the propagating regime, i.e., for 
distances larger than Ath the heat flux is determined by 
Kirchhoff-Planck's law and is limited by the black-body 
value. Note, that the heat flux is dominated by the p- 
polarized modes for distances smaller than 100 nm and 
larger than 10 /im, whereas for distances in between it is 
dominated by the s-polarized modes. 

Now, we introduce the inclusions by using the Maxwell- 
Garnett expression in Eq. (fT3|) and ([II]) . We use the 
same filling factor for both materials, so that we have a 
symmetric situation. In Fig. [5] we show the resulting heat 
flux normalized to the values for the two non-pourous 
SiC plates shown in Fig. HJ We find that for distances 
smaller than 100 nm and larger than 1 /im the heat flux 
becomes larger when we add air inclusions, whereas for 
intermediate distances the heat flux is reduced. 

In order to see how the s- and p-mode contribution is 
changed by the porosity, we show in Fig. [6] (a) and (b) the 
plots for the separate contributions of s- and p-polarized 
modes. It is clear that the p-polarized part of the heat 
flux gets enhanced for all distances when compared to 
the isotropic case, regardless of the filling factor. The 
s-polarized part in turn gives a larger heat flux for dis- 
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Figure 4: Heat flux between two SiC plates over distance with 
Ti = 300 K and T2 = K. The flux is normalized to the value 
for two black bodies Sbb = 459.6 Wm -2 . The contribution 
of the s- and p-polarized part are shown as well. 
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10" 8 10" 7 10" 6 10" 5 10" 4 
d (meters) 

Figure 5: Heat flux between two porous SiC plates over dis- 
tance with Ti = 300 K and T2 = OK. The flux is normalized 
to the value for two SiC plates shown in Fig. [H 

tances larger than about 1 /im and a smaller heat flux 
for distances smaller than 1 /im. Therefore, the smaller 
heat flux found in Fig. [5] for intermediate distances is as- 
sociated to the dominance of s-polarized modes in that 
distance regime. 

In summary, by introducing inclusions we find for large 
and small distances an increase of the heat flux. For the 
propagating regime (d > Ath) this can be understood 
from a simple argument: the vacuum holes simply dilute 
the material so that, according to Kirchhoff's law, the 
reflectivity is decreased and hence the emissivity is in- 
creased. In fact, for / = 1 one would retrieve the black 
body result, since in this case the reflectivity is zero. On 
the other hand, there is no such simple argument for the 
increased heat flux in the near-field region. Here, it is 
necessary to study how the coupled surface modes, which 
give the main contribution to the heat flux for distances 
smaller than 100 nm, are influenced by the introduction 
of the inclusions. This will be done by inspection of the 
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Figure 6: As in Fig [5] but for the (a) s- and (b) p-polarized 
contribution only. 



transmission coefficient in the next section. 



VI. TRANSMISSION COEFFICIENT 

As mentioned before, for the small distance regime 
(d < 100 nm) the heat flux between two isotropic semi- 
infinite SiC-bodies is solely dominated by the p-polarized 
contribution. This remains true for the porous SiC bod- 
ies. In fact, the dominance of the p-polarized contribu- 
tion becomes even greater with increasing filling factors. 
Hence, to understand the observation that by introduc- 
ing some porosity the heat flux becomes larger, it suffices 
to study the p-polarized contribution. 

In Fig. [7] we show the transmission coefficient 
T p (u,K;d) in the (cj, ft)-plane for different filling factors 
and a distance d = 100 nm. In Fig. [71 (a) one can see 
T p (cj, K\d) for two isotropic SiC plates. Here, T p (cj, ft; d) 
is one or close to one for the propagating modes, the total 
internal reflection modes and the coupled surface phonon 
polaritons. In the plotted region one can mainly see the 
coupled surface phonon polaritons, which are responsi- 
ble for the large heat flux at small distances. Now, for 
/ = 0.1 one can see in Fig. [7] (b) that a second coupled 
surface mode appears due to the air inclusions. In addi- 



tion, the coupled surface mode of the bulk SiC is shifted 
to smaller frequencies. When increasing the filling factor 
[Fig. [71(c) and (d)] the upper coupled surface modes shift 
to higher frequencies and become less important for the 
transmission coefficient. On the other hand, the low fre- 
quency surface modes shift further to lower frequencies. 
Between the two coupled surface mode branches a band 
of frustrated internal reflection modes is formed which 
gives also a non-negligible contribution to the transmis- 
sion coefficient. 

In order to get further information we now consider the 
spectral mean Poynting vector (5^) defined in Eq. ([7]) 
for p-polarization only. We have plotted this quantity in 
Fig. [8] at the same distance as before, i.e., d = 100 nm, 
and again for different filling factors. As in Fig. [7] one 
can see the strong contribution of the two coupled sur- 
face mode resonances, which are shifted in frequencies 
when changing the filling factor. Furthermore, one can 
now observe, that when increasing the filling factor the 
low frequency resonance is not only shifted to smaller 
frequencies, but the resonance is also getting stronger. 

The study can now be completed when considering the 
mean transmission factor for the p-polarized modes, that 
was introduced in Ref£ as 



t p(^) = —o I duf(u)T p (u,K-,d) 



(24) 



with u = hw/k^T and / = u 2 e u /(e u — l) 2 . It represents 
the mean transmission coefficient of a mode specified by 
it's wave vector n for a given temperature T and a small 
temperature difference AT between the two bodies. By 
means of this quantity the heat flux can be rewritten in 
a Landauer-like form 9 



(Sz) 



7T 2 fc|T 

3 h 



'^ K T P ( K )AT. 



(25) 



Note, that for nd ^> 1 and n > uo / c the transmis- 
sion coefficient T p (uj,K;d) is exponentially damped [see 
Eq. (|8J)] and therefore also the mean transmission factor 
Tp(ft). This damping determines the wave vector cutoff 
and hence the number of states contributing to the heat 
flux. 

Now, in Fig. [9] we plot T p (k) for a given distance of 
d = 100 nm and different filling factors normalized to 
the mean transmission factor for two semi-infinite SiC 
bodies. For / = 0.1 the mean transmission coefficient 
for the porous SiC increases for intermediate n but de- 
creases for very large k. The increased mean transmis- 
sion factor is due to the second coupled surface mode and 
the frustrated modes, whereas the lower value for large 
wave vectors can be attributed to a stronger cutoff in the 
transmission coefficient, which means that the number 
of contributing modes is decreased. The enhancement of 
the transmission factor due to the surface mode prevails 
and leads to an enhanced heat flux at that distance. The 
same mechanism is responsible for the enhanced heat flux 
for / = 0.3. On the other hand, for larger filling factors 




Figure 7: Transmission coefficient T p (cj, k; d) in the (lj, ft)-plane for two porous SiC slabs with different filling factors (a) / = 0, 
(b) / = 0.1, (c) / = 0.3, and (d) / = 0.5. The distance is fixed at d = 100 nm. 




Figure 8: Spectral mean Poynting vector (S u ) defined in 
Eq. (0 for two porous SiC slabs with different filling factors 
/ = 0,0.1,0.3,0.5 considering only the p-polarized contribu- 
tion. The distance is fixed at d = 100 nm. 



Figure 9: Mean transmission coefficient defined in Eq. (|24[) for 
different filling factors normalized to the isotropic case (/ = 
0). The distance is fixed at d = 100 nm and the temperature 
at T = 300 K. 



the curves change slightly for intermediate k compared 
to the curve for / = 0.3. The contribution in that inter- 
mediate region is due to the second coupled surface mode 
branch and the frustrated modes. But for very large k 
the mean transmission coefficient increases compared to 



/ = 0.3. This means that by introducing a higher poros- 
ity we soften the cutoff of the transmission coefficient. 
Hence, the number of modes contributing to the heat 
flux is increased and results for large filling factors in a 
further enhanced heat flux. 
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The dependence of the cutoff on the filling factor for 
large n can easily be discussed for the transmission coef- 
ficient T p (ft,cj;<f). It was found in Ref£ that the cutoff 
region, i.e., where T p (ft,u;; d) is exponentially damped, is 
given by 



> log 



Im(e) 



1 

2d 



(26) 



when considering two isotropic semi-infinite bodies at the 
surface mode resonance frequency [see also Refs. 36 ]. For 
the uniaxial anisotropic case as considered here, this re- 
lation changes to 



>log 



Im(y^[eX) / 2d 



(27) 



where the permittivities have to be evaluated at the 
surface mode resonance frequency of the semi-infinite 
anisotropic body (see Appendix B). In Fig. [10] we show 
a plot of ftuni/^iso over the filling factor. It is seen that 
by introducing the air inclusions the cutoff first decreases 
and then monotonically increases. This is the same qual- 
itative behavior as observed for the mean transmission 
factor Tp(ft) in Fig. [9] for nd ^> 1. This reasoning con- 
firms that the number of contributing modes is the main 
mechanism for increasing the heat flux at small distances 
and large filling factors (/ > 0.3). 
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Figure 10: Plot of the normalized cutoff value ft'uni/ ^mso over 
filling factor /. 



can be significantly larger than those traditionally mea- 
sured between two isotropic materials in the same non- 
equilibrium thermal conditions. For small filling factors 
we have determined that this enhancement stems from 
additional surface waves arising at the uniaxial material- 
vacuum interface, clearly indicating that such increase 
is intrinsically connected to anisotropy. Indeed, we did 
calculations for isotropically rarifled SiC plates with low 
filling factors (/ < 0.1) and found that the heat transfer 
modification for is much smaller. In contrast, for larger 
filling factors (/ > 0.3) we have shown that, after a thor- 
ough analysis of the transmission factor, the enhance- 
ment in heat transfer arises mainly from the increased 
number of modes contributing to the flux. 
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Appendix A: Green's dyadic in the gap region 

In order to construct the Green's dyadic in the vacuum 
gap we first start with the Green's dyadic in free space. 
If z > z' Weyl's expansion for the Green's dyadic is 37 

G(r,r') = /_^_^__ e hr(z-z') t (A1) 



(2tt) 2 2 7r 



with 7 r = ^/ uo 2 / c 2 — ft 2 , x = (#,2/)* and k = {k x ,k y y. 
The unit dyadic 1 is the unit dyadic in the polarization 
basis and is defined as 



1 = a+ ®a+ + a+ <g>a+. 



(A2) 



The polarisation vectors for s- and p-polarized waves are 
given by 

a+ = - I k x I and a+ = — k y j r I . (A3) 



M o 



VII. CONCLUSION 

We have presented a detailed study of near and far 
field heat transfer between two flat uniaxial media made 
of polar materials (in our case, SiC) in which cylindrical 
inclusions drilled orthogonally to surfaces are uniformally 
distributed. 

After applying the classical stochastic electrodynamic 
theory to anisotropic materials we have shown that, 
for short distances, the heat flux between such media 



By construction both polarization vectors are orthogo- 
nal. For propagating waves they are also normalized. 
The Fourier component G(k]z,z') of the Greens dyadic 
is defined by 



G(r,r') = J-^^G(^z,z')e 



(A4) 



The above expression for the Green's dyadic represents 
the field of a right going wave at z of a source of unit 
strength placed at z' . If a semi-infinite medium is located 



9 



at z > d then this wave will be reflected so that the 
Green's dyadic G(n; z, z') reads at z > z' 

G a (k; z, z') = — \te^ z - z,) + e 2i7rd e i7 ^ +2/) R 2 l (A5) 

27r 



where we have introduced the reflection matrix 



r 2 = e 

M = { s >p} 



r i,j &i 



(A6) 



with the reflection coefficients r 2 j and the polarization 
vectors a~ = —c/(hiuj)(k x j r ,kyj r ^ 2 ) t and k~ = a+. If 
there is now a second semi- infinite medium at z < with 
a reflection operator defined as 



M = {s,p} 



(A7) 



the waves in that cavity will be multiply reflected at the 
boundaries at z = and z = d so that 23 



G a (k;z,z') 



2 7r _ 



le i lr (z-z') + e 2i 7l -d e -i 7r ( z +z') R2 

e 2i 7r <i e i 7r ( Z - Z ')] RjlR2 
e 4i 7r d e -i 7r ( Z +z') R2]Ri]R2 + _ 



(A8) 



Summing up all contributions we get 



G a (k;z,z') 



2 7r 



D 12 e i 7r (z-z') 
D 21 R2e 2i 7r£ / e -i 7r(z+ z') 



(A9) 



where we have introduced 

ID 12 = (1 -Riltae 2 ^)- 1 , 
D 21 = (1 -Mie 2 ^)- 1 . 



(A10) 
(AH) 



The expression in Eq. (|A9|) is not yet the complete 
intracavity Green's dyadic, since we have not considered 
the waves which start from z' as left going waves and 
arrive after being reflected at the boundary at z = at 
z > z' . With the same reasoning as for Qa(k>', z, z') we 
find for this contribution 



&b{k;z,z') 



2 7r 



D 12 Rie i7r(2;+2;/) 
•B 21 R 2 Rie 2i7 ^e i7r(2/ -^ 



(A12) 



Finally, the intracavity Green's dyadic is given by the 
sum of Eq. (IM1) and (IXT2]) yielding 



G, 



2 7r 



D 12 [\e^ z - z,) + Rie i7r(2+2;/) 
f D 21 ^R 2 Rie i7r(2/ - 2) e 2i7r(i (A13) 
+ R 2 e 2i7rd e- i7r(z+2/) 



Appendix B: The resonances in an anisotropic 
material 



The precise location of resonances can be analytically 
determined from expression ([22]) by solving 



eiiej, = 1. 



(Bl) 



Frequencies which satisfy this condition are resonance 
frequencies of medium because they correspond to a flat 
region of the dispersion curve in the (cj, n) plane and 
therefore to strong density of states. Using expressions 
(fT3|) and (fT^|) this equation can be recast into the follow- 
ing form 



with 



bei + ce n + d = 0, 



a = (l-/) 2 , 
6 = e,(l-/)(2/+l), 
c=(l + /)( e 2 /-l), 
d = e,(/-l). 



(B2) 



(B3) 
(B4) 
(B5) 
(B6) 



The solutions of this equation can readily been calcu- 
lated using the Cardano's method^. When inclusions 
are made by pure vacuum (i.e. e\ = 1), these solutions 
are real and read 



1 

- arccos , 
3 V 2 



(- Q 



27TTI 

~3~ 



(B7) 



for n = 0, 1, 2. with 
3ac - b 2 



V : 



3a 2 



27a 2 d-9abc + 2b 3 
and Q = 27^ (B8) 



Only one of these solutions is positive and must be con- 
served to search the resonance frequencies. 
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